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Abstract. The dynamics group of an asynchronous cellular automaton 
(ACA) relates properties of its long term dynamics to the structure of 
Coxeter groups. The key mathematical feature connecting these diverse 
fields is involutions. Group-theoretic results in the latter domain may 
lead to insight about the dynamics in the former, and vice-versa. In 
this article, we highlight some central themes and common structures, 
and discuss novel approaches to some open and open-ended problems. 
We introduce the state automaton of an ACA, and show how the root 
automaton of a Coxeter group is essentially part of the state automaton 
of a related ACA. 
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1 Introduction 

An asynchronous cellular automaton (ACA) is denned in the same manner as a 
classical cellular automaton (CA) in all aspects except the evaluation mechanism. 
As the name suggests, the maps associated to the vertices (or nodes) are applied 
synchronously for a CA, and asynchronously for an ACA. In general, there are 
many ways that one can apply maps asynchronously. For example, one may 
select a vertex at random according to some probability distribution, apply the 
corresponding map, and repeat this procedure. Alternatively, one may select a 
fixed permutation over the vertices and apply the maps in the sequence specified 
by this permutation. This permutation evaluation process would correspond to 
increasing the time by one unit, and would be applied repeatedly to generate 
the system dynamics. An important aspect of having a fixed permutation update 
sequence is that one obtains a dynamical system. This is not necessarily the case 
in the more general situation, such as when the individual states are updated at 
random. 

The analysis of CAs and ACAs does not have the support that the study 
of ODEs has from established fields such as analysis and differential geometry. 
As such, a key goal of CA/ACA research is to make connections to existing 
mathematical theory. We will consider the class of ir -independent ACAs - those 
whose periodic points (as a set) are independent of the permutation update se- 
quence. While this may seem to be a rather exotic property, we have shown that 
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roughly 40% of the elementary CA rules give rise to n- independent AC As [8]. 
Given a ^-independent ACA, one can define its dynamics group. This permuta- 
tion group on the set of periodic points is a quotient of a Coxeter group, and 
it captures the possible long-term dynamics that one can generate by suitable 
choices of update sequence. Its structure can answer questions about the exis- 
tence and non-existence of periodic orbits of given sizes. 

In this paper, we will revisit the notions of Coxeter systems and sequential 
dynamical systems (SDSs). An SDS is a generalization of an ACA (assuming 
a fixed update sequence) where the underlying graph is arbitrary, and is not 
limited to being a regular lattice or circle (i.e., a one-dimensional torus). We will 
show how the word problem for Coxeter groups is related to functional equiva- 
lence of SDS maps. This forms the basis for our next result, on how conjugation 
of Coxeter elements corresponds to cycle equivalence of SDS maps, and addition- 
ally, how this extends from conjugacy classes to spectral classes. After defining 
dynamics groups and showing how they arise as quotients of Coxeter groups, 
we show how key features of mathematical objects in both the fields of SDSs 
and Coxeter groups are encoded by finite (or infinite) state automata. We illus- 
trate this by explicit examples, and then close with a table summarizing these 
connections. 

2 Background 

A Coxeter system is a pair (W, S) consisting of a group W generated by a set 
S = {si, . . . , s n } of involutions given by the following presentation 

W=( Sl ,...,s n \sf = l, (s iSj ) m ^ = l), 

where m(si,Sj) > 2 for i ^ j. Let S* be the free monoid over S, and for each 
integer m > and distinct generators s,t G S, define 

(s, t) m = stst ■ ■ - t G S* . 

m 

The relation (s, t) m ( Si t) — (t, s) m ( st ) is called a braid relation, and is additionally 
called a short braid relation if m(s, i) = 2. Note that s and t commute if and 
only if m(s,t) = 2. A Coxeter system can be described uniquely by its Coxeter 
graph T, which has vertex set V — {1, . . . ,n} and an edge {i, j} for each non- 
commuting pair of generators {sj, Sj}, with edge label m(sj, Sj). 

Switching to ACAs and SDSs, let r be an undirected graph (called the 
base graph or dependency graph) with vertex set V = {l,...,n}. We equip 
each vertex i with a state x% G K where AT is a set called the state space, and 
a vertex function fi that maps (or updates) Xi(t) to Xi(t + 1) based on the 
states of its neighbors (itself included). Unless explicitly stated otherwise, we 
will assume that A = F2 = {0,1}, which is the most commonly used state space 
in cellular automata research. If the vertex functions are applied asynchronously, 
it is convenient to encode /» as a F -local junction Fi : K n — > K n defined by 

Fi{x\, . . . , X n ) — {x\ , . . . , Xi—\ , fi(x\ , . . . , X n ), X^-j-i , . . . , X n ) . 
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If fi does not depend on all n states, it may be convenient to omit the fictitious 
variables. Given a sequence of local functions and a word w = W1W2 ■ ■ ■ w m G V* 
called the update sequence, the SDS map F w is the composition of the local 
functions in the order prescribed by w, i.e., 



F„ : K r 



K' 



F = F o F 



■■°F W2 o F W1 



SDSs represent a generalization of ACAs, which are usually defined over a regular 
grid, such as Z or Z„. The following example illustrates some SDS concepts - 
see [14] for a more complete treatment. 

Example 1. We take r — Circ4 as base graph (see Figure 1) and use K — 
{0,1} as the state space. Also, we take all vertex functions to be Boolean 
nor-functions given by nor: K 3 — > K where nor(cc, y, z) equals 1 if x = y = 
z = and otherwise. In this case we have, for example, Fi(x±, X2, 2:3, £4) = 
(nor(iE4, sci, X2), X2, £3, Xi). Using the update sequence ir = 1234, we obtain 

^(0,0,0,0) = (1,0,0,0) 
F 2 oFi(0,0,0 ! 0) = (l,0,0,0) 
F 3 oF 2 oF 1 (Q, 0,0,0) = (1,0, 1,0) 
Fi 0F3 oF 2 oFi(0,0,0,0) = (1,0,1,0) , 

and thus F n (Q, 0, 0, 0) = (1, 0, 1, 0). The phase space of the map F n is the directed 
graph containing all global state transitions and is shown in Figure 1. 
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Fig. 1. The base graph _T = Circ4 and the phase space of F n from Example 1. 



3 The word problem 

A fundamental question given any finitely presented group (S \ R), is when do 
two words 

w = W1W2 ■ ■ ■ w m , and w' = w^w^ ■ ■ ■ w k 

in S* yield the same group element? This is the word problem, and it is in gen- 
eral undecidable. However, there are many classes of groups for which the word 
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problem is solvable. A classic result in Coxeter groups, known as Matsumoto 's 
theorem [5, Theorem 1.2.2], says that any two reduced expressions for the same 
clement differ by braid relations. Matsumoto's theorem provides an algorithmic 
solution to the word problem for Coxeter groups. 

There is an analog of the word problem for SDSs. Specifically, given two 
update sequences w, w' £ V* , when are the corresponding SDS maps 

F w = F Wm o F Wm _ 1 o • • • o F W2 o F W1 , F w , = F w > k o F w > k _ i o • • • o F w > 2 o F w[ 

equal as functions, or equivalently, when do they have identical phase spaces? 
This is clearly solvable because there are only finitely many functions F£ — > F£ . 
However, it would be desirable to solve this problem algorithmically for general 
SDSs, without resorting to checking the image of all 2™ global states. 

4 Equivalences on Dynamics and Acyclic Orientations 

In this section, we show how topological conjugation of SDS maps corresponds 
to conjugation of elements in a Coxeter group, and how this connection leads 
to a coarser equivalence relation when the graph r has non-trivial symmetries. 
Acyclic orientations are mathematically convenient to capture several types of 
equivalences on permutation SDS maps, as well as on Coxeter elements in Cox- 
eter groups. A Coxeter element is the product of the generators of S in some 
order. Every Coxeter element defines a partial ordering on S, which we can rep- 
resent by an acyclic orientation of r. Specifically, for a Coxeter element c, define 
the orientation (T, c) so that edge {i,j} is oriented (i,j) if s» appears before Sj 
in c. It is easy to show that this is well-defined, and that it defines a bijection be- 
tween the set Acyc(-T) of acyclic orientations of r and the set C(W) of Coxeter 
elements of W . 

Next, consider conjugating a Coxeter element c = s Xl ■ ■ ■ s Xn by the initial 
letter s — s Xl , which results in a cyclic shift of the word: 

SCS — S Xl (s Xl S X2 • • • S Xn *)S Xl — $X2 S X 3 ' ' ' $x n &Xi • 

The corresponding acyclic orientations (r, c) and (r, scs) differ by converting 
the source vertex of (_T, c) into a sink. This source-to-sink conversion generates 
an equivalence relation ^ K on Acyc(_T), and it was recently proven (see [4]) that 
(r, c) ^ K (r, c') if and only if c and d are conjugate. (Note that the "if" direction 
is obvious; the "only if" direction is difficult). 

Turning to SDSs, let S n C V* be the set of words where each vertex appears 
precisely once, which we may identify with the permutations of V. Each permu- 
tation 7r £ S n defines a partial ordering on V, and there is a natural map from 
Acyc(/ n ) to the set of permutation SDS maps (it is mapped to F v ). Two finite 
dynamical systems <fi, ip : K n — > K n are said to be cycle equivalent if for some 
bijection h: K n — > K n we have V , |per(v) ° ^ — ^°0|p C r(0) ; where Per(0) denotes 
the set of periodic states of <j>. The following result provides the connection be- 
tween K-equivalence of acyclic orientations and cycle equivalence of permutation 
SDS maps. 



Coxeter Groups and ACAs 5 



Theorem 1 ([12]). If (r,n) ^ K (r,a), then F„ and F a are cycle equivalent. 

If the automorphism group Aut(-T) is non-trivial, we can say even more. The 
group Aut(T) acts on Aut(.T)/~ K by 7 • {(r, 7r)] = {(r, 771")], which gives rise 
to the equivalence relation ~ R on Aut(-T)/^ K . This coarser equivalence relation 
also has an interpretation in the settings of both Coxeter groups and SDSs. 

If (W, S) is a Coxeter system with IS"! = n, let V be an n-dimensional real 
vector space with basis {cci, . . . , ot n }. Put a symmetric bilinear form B on V, 
defined by B{an, ctj) = — cos(n/m(si, Sj)) . The group W acts on V by 



and the set of elements <P — {won \ w £ W ,i = 1, . . . ,n} are called roots. This 
action is faithful and preserves the bilinear form B. Geometrically, the root SiV 
is the reflection of v across the hyperplane aj~, and so there is a representation 
p: W — > GL(V), defined on the generators by 



called the standard geometric representation of W (see [1,7]). This allows us to 
view elements in W as matrices, and hence we can speak of the characteristic 
polynomial of any given w £ W. 

Now, if (r, c) and (-T, c') differ by some 7 G Aut(T), then p{c) and p(c') 
are similar as linear transformations. Specifically, they are conjugate in GL(V) 
by the permutation matrix P 1 of 7. In this case, we say that c and d have the 
same spectral class, because p{c) and p(c') have the same multiset of eigenvalues. 
Clearly, this is a weaker condition than conjugacy, and so all Coxeter elements 
in the same K-equi valence class have the same spectral class. 

Similarly, in the context of SDSs, if (r, ir) ^ R (T, a), then the SDS maps 
and F (j are cycle equivalent, due to the following argument. If 7 £ Aut(_T), then 
the permutations 7r and 77r give topologically conjugate SDS maps, F^ and F 17T . 
Strictly speaking, this requires the maps f t to be Aut(-T)-invariant (see [12]), 
a condition which is frequently satisfied in practice, such as when all vertices 
of the same degree share the same symmetric function (e.g., logical AND, OR, 
Majority, Parity, threshold functions, etc.). Since topologically conjugate maps 
are cycle equivalent, our statement follows. 

It is worth mentioning the role of the Tutte polynomial here [15]. The Tutte 
polynomial of a graph r is a 2-variable polynomial Tr(x,y) that satisfies a re- 
currence under edge deletion and contraction, and plays a central role in graph 
theory. Many graph counting problems are simply the evaluation of the Tutte 
polynomial at some (xo,yo) £ZxZ. For example, | Acyc(F)| = Tr (2,0) and 
I Acyc(_T)/^ K | = Tp(l, 0). Thus, Tp(2, 0) counts the number of Coxeter elements 
in the Coxeter group with Coxeter graph r, and it bounds the number of per- 
mutation SDS maps in an SDS with dependency graph r. This bound is known 
to be sharp for certain classes of functions [14]. Similarly, Tp(l,0) counts the 
number of conjugacy classes of Coxeter elements (see [4, 10]), and it bounds the 
number of cycle equivalence classes of SDS maps (see [12]). 



Si : v M> v — 2B(v, ati)ati 



(1) 




(2) 
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5 Groups 

A sequence F = {Fi, ■ ■ ■ ,F n ) of local functions is tt -independent if Pei(F 7T ) = 
Per(F CT ) for all ir, a e S n . Note that this is an equality of sets; we do not assume 
anything about the organization of the respective periodic points into periodic 
orbits. In this case, each Fi permutes the periodic points, and these permutations 
generate the dynamics group of F, denoted T>Q(F). Let F* denote the restriction 
of ^ to Per(F T ). Because Fi only changes the i th coordinate of a state, and since 
we assume that K = ¥2, F* oF* is the identity function on Per(i 7 ' T ). If we define 
rriij :— \F* o Fj |, then there is a surjection 

(s u ...,s n I sj = 1, {s iSj ) m ^ = 1) — > VG(F) , (3) 

showing that dynamics groups are quotients of Coxeter groups. The particular 
homomorphism is determined by adding relations to the presentation of the 
Coxeter group, and these relations arise because the state space is F2. Thus, 
dynamics groups are in a sense "reflection groups over F2." An open-ended 
research problem is to give an efficient presentation of the dynamics groups of 
an SDS based on the functions, i.e., to determine these extra relations. 

When the base graph r of an SDS is the circular graph Z„, and the local 
functions are all identical, the resulting SDS is an elementary AC A. Each local 
function Fi takes {ajj-i, Xi, Xi+i} as input, and is completely described by the 
following rule table 





111 


110 


101 


100 


Oil 


010 


001 


000 


fi {p^i— 1 ) ^ii ^i+l) 


a 7 


o 6 


a 5 


(24 


«3 




ai 


a 



Clearly, there are 2 2 ' = 256 such choices of functions, which can be indexed 
by k = ^2 a i^ % e {0: • • • j 255}. The corresponding sequence of local functions is 
denoted ECAfe. In [8], it was shown that ECA^ is n- independent for precisely 
104 values of k. Moreover, this holds for all n > 3. The dynamics groups of these 
104 rules were classified in [11]. Among some of the interesting groups were 
£>g(ECA 60 ) = SL„(F 2 ) and Vg{ECA k ) = Z% for k £ {28,29,51}. Moreover, 
other dynamics groups were found computationally to be either the symmetric 
or alternating groups, with the size depending on the n th Fibonacci or Lucas 
number, leading to a few conjectures. 

6 The root automaton 

The dynamics of all possible SDSs given a sequence of F-local functions F = 
(Fi, . . . ,F n ) can be encoded by the state automaton of the sequence. This is 
a directed graph <P with vertex set K n - the set of global system states, and 
directed edges (x,Fi(x)) for each x £ K n and each i G V. Label such an edge 
with the index i corresponding to its vertex function; see Figure 2 for an example. 
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Fig. 2. The state automaton # for the SDS in Example 1. Horizontal/vertical dashed 
lines and arrows indicate horizontal/vertical wrap-around, and arrowheads are omitted 
from the bidirectional edges for clarity. 



The image of a state x G K n under an SDS map F„, where tt = ~k\~ki ■ ■ ■ 7r„, 
is represented on the state automaton by a path in <P. Specifically, start at vertex 
x and traverse the path 

1 -^7Ti (■£) ; ^7T2 ^7Tl (*^) 5 ' * * 7 ^7T n ' ' ' -^7T2 -^7Ti (^) — -^7T (^) • 

The phase space of F„ can be easily derived from the state automaton - it is 
the graph with vertex set K n and an edge (x, y) for every directed path from 
a state x to y that traverses a path of edges labeled 7ri,7T2, . . . ,7r n . Note that if 
F is 7r-independent, then T>Q(F) acts on Per(F). In this case, all of the edges 
within Pei(F) are bidirectional, and so we may view them as undirected. 

This is the SDS analog of the action of W on <P C V, as described in (1). Since 
V is any n-dimensional vector space, we can identify it with R™, and assume that 
the basis elements are ex, = e,, the standard unit normal vectors. This associates 
roots with vectors in R™, and we partially order <P by < componentwise (z < z 1 
iff Zi < z[ for each i) to get the root poset. It is well-known that for every root, 
all non-zero entries have the same sign, thus we have a notion of positive and 
negative roots, and the root poset has a positive side <P + and a negative side, 
<P~ , with <P = <P + U <P~ . The image of Si under the geometric representation 
from (2) is a linear map Fi : R n — > R™, where 

n 

Fi - {z\,-.-,z n ) i — > (z lt . . . , Zj + ^ 2 cos{-K/m itj )zj, z i+ i, ...,z n ) . (4) 

i=i 

To summarize, Fi changes the i th entry of a vector by flipping its sign and then 
adding each neighboring state Zj weighted by 2 cos^/my). 

In 1993, Brink and Howlett proved that Coxeter groups are automatic [2], and 
soon after, H. Eriksson developed the root automaton [3]. The root automaton 
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has vertex set <P and edge set {(z, s^z) z e ^, Sj e 5}. For convenience, label 
each edge (z,Sjz) with the corresponding generator Sj. It is clear that upon 
disregarding loops and edge orientations (all edges are bidirectional anyways), 
we are left with the Hasse diagram of the root poset. We represent a word 
w = s Xl s X2 ■ ■ ■ s Xm in the root automaton by starting at the unit vector e Xl € <P + 
and traversing the edges labeled s X2 , s X3 , . . . , s Xm in sequence. Denote the root 
reached in the root poset upon performing these steps by r(W, w). The sequence 

e Xl = r(W, s Xl ) , r(W,s Xl s X2 ), ... , r(W, s Xl s X2 ■ ■ ■ s Xm ) = r(W,w) , 

is called the root sequence of w. If r(W, s Xl s X2 ■ ■ ■ s Xi ) is the first negative root 
in the root sequence for w, then a shorter expression for w can be obtained by 
removing s Xl and s Xi . By the exchange property of Coxeter groups (see [1,7]), 
every word ai£S* can be made into a reduced expression by iteratively removing 
pairs of letters in this manner. Thus, the root automaton can algorithmically 
detect reduced words. 

We conclude with an example that illustrates these concepts, and shows how 
the root automaton of a Coxeter group is essentially a connected component of 
the state automaton of an sequential dynamical system with state space K = R. 

Example 2. Let W = if 4, which has Coxeter graph as shown in Figure 3, and 
presentation (using a, b, c, d instead of s 1; s 2 , s 3 , s 4 ): 

Hi = (a,b,c,d I a 2 , b 2 , c 2 , d 2 , (ab) 5 , (be) 3 , (cdf , (ac) 2 , [ad) 2 , (bdf) . 

It is well-known (see [7]) that if 4 is a finite group of order 14400, and is the 

5 



abed 

Fig. 3. The Coxeter graph r of the group W = H 4 . As is customary, edge labels of 3 
are suppressed. 

isometry group of the 120-cell and its dual, the 600-cell, two of the six regular 
4-polytopes. Thus, the root poset consists of 14400 roots. A portion of the 
root automaton is shown in Figure 4. Recall that the root automaton is built on 
top of the root poset - stripping away the self-loops and edge labels leaves the 
Hasse diagram of <P. The dotted-line in Figure 4 shows the boundary between 
the positive roots <P + and negative roots <&~ . The non-loop edges of the root 
automaton are all bidirectional - arrowheads are omitted for clarity. 

Consider the word w = abdeabaebca e H 4 . Starting at e a = (1,0,0,0) (see 
Figure 4), and traversing the edges labeled b, d, c, a, 6, a, c, b, c, a in sequence, we 
see that the first negative root in the root sequence of w is r (W, abdeabaebe) . 
Therefore, removing the first instance of a and the last instance of c from w 
results in bdeabaeba, a shorter expression for w. It is easily checked that no matter 
where we begin in bdeabaeba, the corresponding path in the root automaton 
consists of only positive roots. Therefore, bdeabaeba is a reduced word in if 4 . 
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Fig. 4. Part of the root automaton of the group W = -H4. Here, <j> = 2cos(-7r/5), the 
golden ratio. The dotted line separates the positive roots from the negative roots. 



7 Summary 

This paper presented a collection of results connecting properties of Coxeter 
groups and properties of the dynamics of ACAs/SDSs (see Table 1). These newly 
established connections provide possible avenues for ACA research. In a larger 
setting, we hope that our example linking properties of asynchronous, finite 
dynamical systems and group theory can provide inspiration for other approaches 
seeking to better understand the dynamics of ACAs through the use of existing 
mathematical structures and theory. 
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Coxeter groups 


Sequential dynamical systems 


Graph r <- 


-¥ Coxeter graph 


Dependency graph 


Acyc(r) ^ 
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W = S 7r ( 1 )S 7r (2) • • • S,r( n ) 


Permutation SDS maps 
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K-equiv. <— 
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